Coherent emission from disordered arrays of driven Josephson vortices 
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We propose a mechanism of coherent emission from driven vortices in stacked intrinsic Josephson 
junctions. In contrast to super-radiance, which occurs only for highly ordered vortex lattices, we 
predict resonant radiation emission from weakly correlated vortex arrays. Our analytical results 
for the THz wave intensity, resonance frequencies, and the dependence of THz emission power on 
dissipation are in good agreement with the ones obtained by recent simulations. 
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I. INTRODUCTION 

It has been experimentally observed 1 - - — and con- 
firmed both analytically 4 ^ and numerically 6 : 7 that mov- 
ing Josephson vortices (JVs) emit sub-THz electromag- 
netic radiation. Tcra-Hertz radiation has applications in 
physics, astronomy, chemistry, biology, and medicine 9 . 
This motivates recent proposals^ for THz filter o 10 ' 12 , 
detectors 1 - 3 -, quantum device a 14 ' 15 , and emitters^ based 
on highly anisotropic layered superconductors (e.g., 
Bi2Sr2CaCu20s+5), that can be modelled as coupled in- 
trinsic Josephson junctions (IJJ). 

The ultimate challenge in this field is to produce co- 
herent THz radiation. It is commonly believed that 
this goal can be achieved by controlling super-radiance 
from highly ordered vortex lattice s 16 : 17 . A vortex lat- 
tice is deemed necessary because the constructive inter- 
ference of Josephson plasma waves from individual JVs is 
strongly suppressed by small amounts of disorder. Unfor- 
tunately, driven periodic lattices are often very unstable 
(especially in the presence of impurities, defects, and pin- 
ning centers), and moving JVs form either a mixture of 
coexisting different lattices 1 -^ or even disordered arrays 6 -. 
Moreover, a broad radiation spectrum by individual vor- 
tices results in a broad spectrum of the emitted radia- 
tion (e.g., Ref. Hf), in contrast to a desirable resonant 
IJJ, where coherent radiation is characterized by sharp 
spectral lines. 

In this context, recent interesting simulations by 
Tachiki and coworkers 6 show that coherent radiation may 
be generated by JVs moving as disordered arrays, instead 
of just ordered ones. This raises the question as under 
what conditions JVs in layered superconductors emit co- 
herent radiation. Solving this problem is crucial to the 
effective design of I^S^CaC^Os+a-based THz emit- 
ters. 

The experimental demonstration 3 of THz radiation in 
zero magnetic field and various failed attempts at de- 
tecting THz emission in the presence of magnetic fields 
cast serious doubts on the initial idea that moving JVs 



can radiate in this frequency domain. Indeed, the now 
prevailing interpretation is that JVs ought to be con- 
sidered as perturbing degrees of freedom, which destroy 
the layer coherence and thus cause the suppression of 
THz radiation. In this study, however, we reach the con- 
clusion that, under appropriate conditions, applied mag- 
netic fields do help amplify and tune THz emission. This 
interesting result is also consistent with the recent sys- 
tematic studies in Ref. Ojl 

Below, we show that the nonlocal nature of JVs in lay- 
ered superconductors is responsible for a two-scale dy- 
namics. A longer scale, Xem, characterizes the inter- 
vortex magnetic interaction; spatial dispersion, or dis- 
order, of vortices up to such a scale has no appreciable 
impact on the radiation mechanism. In other words, in 
contrast with super-radiance, which is suppressed by vor- 
tex disorder, in our approach radiation coherence is pre- 
served even if the vortex distribution can become appre- 
ciably modulated by the radiation itself for wavelengths 
shorter than Xem- The shorter relevant length scale 
Xq determines the cross-section of the nonlinear vortex 
core (where the linear approximation sin ip rs ip with the 
gauge invariant phase difference across the junction hav- 
ing a vortex is not valid). According to this picture, 
for Xg <C Xem and for a sufficiently high vortex den- 
sity, radiation is emitted through a linear mechanism, as 
from a JV lattice, whereas the vortex-radiation coupling 
occurs mainly in the inner JV cores. Under these con- 
ditions, the magnetic interaction among vortices is much 
weaker than their interaction with the emitted radiation. 
Our approach explains the spatial modulation of the JV 
density numerically found in Ref. [6||7|. Our analytical es- 
timates, based on a one-dimensional sine-Gordon (sG) 
model, prove to be in good agreement with their simula- 
tions and explain their results. 

Let us now summarize a central idea of our approach. 
Consider a moving JV lattice emitting radiation. This 
radiation will bounce back and forth the sample edges, 
like in a laser cavity. This radiation accumulates and 
creates a standing wave with a wavelength about Xem- 
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This standing wave modulates the JV density which is 
now in resonance with the standing wave. This positive 
feedback enhances the radiation of vortices. Namely, the 
JV motion emits radiation, which is weaker at first. This 
radiation bounced inside the sample (acting as a cavity) 
locks the collective motion of the JVs. This collective mo- 
tion produces stronger emission. The JVs then interact 
more strongly with the electromagnetic standing wave, 
compared with the now much weaker vortex-vortex in- 
teraction. Thus, the triangular vortex lattice, produced 
by the vortex-vortex interaction, is finally replaced by a 
more disordered, but still modulated by the radiation, 
vortex structure. 



II. NONLOCAL SINE-GORDON MODEL 

Layered superconductors can be considered as stacks 
of strongly interacting IJ Js. As the superconducting lay- 
ers are only a few nanometers thick, i.e., the inter-layer 
distance s is much smaller than the magnetic field pene- 
tration depth X a b, the currents flowing through different 
junctions are coupled. On neglecting, for the time be- 
ing, external drives and internal dissipation, a system of 
stacked IJJs is well described by the coupled sine-Gordon 



equation; 
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(n) ^ 

^V+sin^") 



A 2 ^=0, (1) 



where <pW is the gauge invariant phase difference across 
the nth junction. Here, oj p is the Josephson plasma fre- 
quency, A the London penetration depth (A/A a b = 7) 
along the layers, and the operator A 2 is defined by 
A 2 / = /(™+ x ) - 2/W — 

A full analysis of this set of equations is a complicated 
problem which requires numerical simulation. However, 
if we restrict ourselves to the case of moderate mag- 
netic fields, when JV cores do not overlap, we can re- 
duce Eq. (TT|) to an effective ID problem. Indeed, as was 
shown inii (see Eqs. (21) and (22) and Fig. 2 there), the 
phase difference ip decreases very fast away from the junc- 
tion where a vortex located. Thus, a reasonable strategy 
could consist in neglecting the nonlinear couplings be- 
tween junctions at a distance of some s from the vortex 
center; on solving the linearized equations (fTJ) for such 
junctions, one would end up with a few coupled nonlinear 
equations for a few junctions in the vicinity of the vortex 
center. The case when the nonlinearity was restricted 
to one junction only has been considered in Ref. |5f The 
coupled junction system of Eq. (fT]) boils down to a ID 
Josephson junction described by a nonlocal sine-Gordon 
equation. We assume below that retaining the nonlinear 
coupling between more junctions can lead to the same 
nonlocal ID sine-Gordon equation with additional noise- 
like weak perturbations. 

For simplicity, let us consider the pair of adjacent junc- 
tions j and (j : + 1) locating a moving JV. We then reduce 



the description of the IJJ stack to a ID problem by as- 
suming the nonlinear coupling to be important only for 
the paired junctions and linearizing Eq. (TT]) for all other 
junctions (i.e., for n 7^ j, j + 1). It is interesting to note 
that the importance of the interaction between two neigh- 
boring junctions is numerically well established 21 . More- 
over, Koshelev^ 2 - recently reduced the multi-junction sys- 
tem to two coupled junctions, and this model reproduced 
the simulation data in Ref. [2l| and interpreted the exper- 
imental results of Ref. H. 

Following the approach in Refs. I5I I23II24I the equation 
for the averaged phase difference across a junction pair 
(p = (ip^ +l ^ + ip^)/2 can be written as 

r "" '-sin^= ^ / dx'K Vxx { x ') + P[^] sin^, 
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where Kq is the modified Bessel function and 

P[ip] = 1 — cos ip 
with ip = (<^' +1 ) - <pW>)/2. The length 



(2) 



(3) 



defines the size of the JV core. Again, the contribution 
of the next-to-neighbor junctions to the dynamics of the 
tagged JV weakens fastii with their distance from the 
vortex center, thus, allowing all other nonlinear equations 
(fT]) to be replaced by an effective nonlinear medium. 

An additional equation for ip can be derived for a pair 
of JVs in two adjacent junctions, so that the equations 
for ip and ip form a closed set^^. However, when extend- 
ing Eq. ((2|) to describe the collective motion of N travel- 
ling JVs (randomly distributed along Ni stacked IJJs of 
length L), the phase <p can be regarded as a mean-field 
superposition of the n = N/Ni individual JVs phases, 

ipv , contained in one layer, only. 

In our one-IJJ description we assume that ip is rela- 
tively small; this may be the case, for instance, due to 
the random superposition of the vortex dynamics in dif- 
ferent junctions. Anyway, the good agreement between 
the analytical results reported here and earlier numer- 
ical simulations, validates a posteriori our assumption. 
Thus, the functional P[ip] can be modelled as a spatial 
perturbation S + eP(x), where the real function P(x) can 
be either periodic or random in x, depending on the op- 
erating conditions. The constant e is a measure of the 
strength of the perturbation, while the small offset 8 can 
be conveniently eliminated by rescaling the dimensional 
parameters uj p and Ag, as appropriate. Here, we assimi- 
late such a perturbation as an effective quenched Gaus- 
sian disorder along the IJJs; that is, P{x) is modeled as 
a random, delta-correlated function with 



(P(x)) = 0, (P(x)P{x')) = 2S{x - x'), 



(4) 



and (...) denoting the average over different disorder re- 
alizations. The constant e will be taken as a perturba- 
tion parameter and only effects to leading order in e will 
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be considered. Moreover, deviations from the Gaussian 
statistics, implicit in the definition of P[^>], are assumed 
to be negligible within this approximation. 



A. Josephson vortex array 

We now introduce dimensionless units by expressing x 
and t in units of the characteristic length Xem and the 
reciprocal of the plasma frequency cj p , respectively; that 
is 

x — )■ x — x/Xemj £—>•£ = u p t. 

As a consequence, the system characteristic lengths A and 
Xg get rescaled as follows: 

X — > X = X/Xem, Xq — > Xq = Xem/X = 1/ A, 

and, of course, Xem ~ > Xem = 1- Correspondingly, 
tp{x, t) is given in units of the magnetic flux quantum 
4>o and all speeds in units of oj p X. Hereafter, for the sake 
of simplicity, we shall only use dimensionless variables 
and, therefore, omit the "tilde" notation altogether. 

The field tp(x, t), corresponding to a dense distribution 
of JVs traveling with speed V <C 1, can be expanded 
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<p(x, t) = p (x — Vt) — k sm[p (x — Vt)] + 



(5) 



where p — 2irp, and p = n/ L denotes the linear JV den- 
sity with number n of vortices located along the length L. 
The linear term in Eq. ([5]) corresponds to the phase dif- 
ference of a uniform vortex spatial distribution, whereas 
the periodic correction accounts for a residual phase mod- 
ulation on the lattice scale 1/p, with amplitude k to be 
determined self-consistently. In the expansion ([5]) we as- 
sume high JV densities, p » 1, and small amplitudes 
k. On inserting expansion ([5]) for (p(x, t), the nonlocal 
field equation ([2} can be approximated to an effective 
sine-Gordon equation, where 

If ( \x — x'\ \ . . 2 

— / K I — I tp xx (x ) dx -> c p tp xx , (6) 



and, consistently, 



with 



c% = [1 + (pA) 2 



K = ijvlCppf 



In the regime considered in Refs. H@, where A ^> 1 and 
7v — 1, the parameter c p can be further approximated 
to (pX)~^. 

The validity condition for truncating the expansion ([5]) 
to its first order, k <C 1, or equivalently c p p ^> 1, implies 



a direct core-core interaction; that is 1/A >• 1/p. We 
recall that here 1/A represents the size of a vortex core 
in dimensionless units. Accordingly, for A > 1, c p must 
be regarded as the maximum velocity of a vortex array 
in a layered superconductor, to be compared with the 
maximum dimensionless velocity 1/A of a single vortex 
(i.e., ui p Xg in dimensional units 

In the opposite limit, pjX <C 1, vortices only weakly 
interact on the magnetic length scale Xem (rescaled here 
to 1); the limiting velocity c p grows larger than 1/A and 
the JV array becomes unstable. 



B. Radiation mechanism 

The emission of radiation by fast moving JVs also takes 
place on the magnetic length scale Xem- The effective 
phase difference ip associated with an array of JVs moving 
along an IJJ, thus, obeys the perturbed local sine-Gordon 
equation 

tptt - c p 2 tp xx + sin ip = —P<p t - f + eP(x) sin (p. (7) 

Note that, in leading order, tp = p (x—Vt), as can be seen 
from Eq. ([5]). Here, for completeness, we have restored 
the viscous term — 0<pt and the current-induced drive /, 
that allow us to control the net JV speed V (see, e.g., 
Ref.H). 

Like in the more conventional single sine-Gordon- 
soliton perturbation scheme o 27 i 28 , we consider the Ansatz 



<p{x,t) -> tp(x,t) + x(x,t), 
which, inserted in Eq. ([7]), yields 26 



(8) 



Xtt - c p 2 Xxx + (cos ip)x = - fixt + eP(x) sin ip, (9) 

where tp = is the ground state and only terms 0(e) have 
been kept. The wavenumber q and the angular frequency 
uj of the unperturbed plasmon modes (i.e., for (3 = e = 0) 
form a continuum spectrum 2 ^, with 



2 1, 2 2 

L) = 1 + c p q . 



(10) 



However, as for the field (|5|) with p » 1, the radiation- 
vortex coupling (costp)x becomes negligible, and the 
plasma wave dispersion relation can be approximated to 
u) = c p \q\. On introducing the spatial Fourier compo- 
nents of P{x) and x(a;,i), defined by 



P(x) 



P(k) 



e tkx dk. 



xOM) = 



xSY qX dq , 



Eq. © can be rewritten a; 



.2!) 



dt 
where 



e 

2i. 



B(q) = (x q -i\q\ Xq y 
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is directly related to the spectral density of the array 
emission power, 

7r at 



that is 



W(q) 



2c 



2 r 



\P(q-p)\ 2 (13/2) 
((3/2) 2 + (w-pVy 



+ {p^ -P) 



(11) 



Here we use that P(x) is a real function, so that P*(k) = 
P(—k). The notation "(p — > —p)" denotes the symmetric 
term obtained by replacing p — > — p in the first term 
inside the square brackets. This means that two waves 
propagate in opposite directions with the same frequency 
u\ for \P(q — p)\ = \P(q + p)\, they generate standing 
plasma oscillations, like those reported in Ref. [HQ- To 
simplify our notation, hereafter we restrict ourselves to 
JVs driven in one assigned direction, say V > 0. 

The spectral emission power (s.e.p.) (fTTT) is key to our 
analysis of a resonant IJJ. The spectrum W(q) can be 
easily specialized for any choice of P(x). In the case of 
quenched Gaussian disorder, see Eq. (HJ, 



(\P(k)\ 2 ) = -, 



(12) 



so that on disorder-averaging Eq. (fTTj) we obtain the IJJ 
spectral emission power per unit of length 



e 

4tt 



(3/2 



(p/2y + {u~pvy 



(13) 



This spectrum holds for j3 <C oj or, equivalently, for V S> 
(3/p, and has a sharp resonance maximum 



for 



e 

2^/3 



(14) 



(15) 



C. Vortex dynamics 

Subject to a drive / produced by an externally-applied 
electrical current, the vortices in an IJJ flow with an av- 
erage speed V and, simultaneously, their cores interact 
with the electromagnetic waves they radiate. For the rel- 
atively weak vortex core repulsion, 1/A > 1/p, simulated 
in00, we expect that the vortex array can be modulated, 
both in space and time, by the resonant plasma modes. 

To express the average speed V of a JV array with 
p/A > 1 as a function of the drive /, from Eq. (J7J we 
derive the energy balance equation per unit of length of 
radiating IJJ2S 



Equation (| 16[) tells us that the rate at which the drive 
pumps energy into the system (right-hand-side), must be 
equilibrated by the radiative, w(V), and the viscous loss, 
(3(pV) 2 , of the soliton array <p(x, t) (left-hand-side). 

For V ^> (3/p, the total emission power of the radiating 
sine-Gordon solitons, 



2 ' 



(17) 



is computed by integrating the spectral emission power 
(fTTj) ; solving the ensuing Eq. (fT5|) with respect to V, we 
obtain 



V(f) 



/±(/ 2 -2/3e 2 )^ 
2(3p 



(18) 



where only the rising branch with the + sign is stable 2 ^. 
Therefore, the observable velocity- drive characteristic 
V(f) is expected to show a step at 



/th = (20)*e 



(19) 



and to grow linearly with / for / 3> /th, when the radi- 
ation loss becomes negligible, namely 



V = —. 

f3p 



(20) 



Note that, at variance with an emitting JV lattice 2 ^, no 
multiple hysteretic steps in the V(f) are predicted. In- 
deed, the condition / > f t h simply implies that the ef- 
fective phase ip is not pinned by disorder 32 ; for f < fth, 
instead, the JV array can move only by creeping, namely, 
through the nucleation and the subsequent migration of 
array defects^i^. Creeping is likely responsible for the 
smooth low-current J-V characteristics shown in Fig. 4 
of Ref. 6. Moreover, in the linear regime (|2T)j) the wave- 
lengths A r of the emitted radiation are expected to be 
much shorter than the length L of the IJJ (sec below), 
so that corrections due to the appropriate standing-wave 
periodic boundary conditions are of the order of \ r /L. 

A vortex is sensitive to the radiation field only when 
the wavelengths A,- excited in the IJJ are larger than its 
size. For the parameter choice of Refs. 00, this can only 
occur on the JV core scale Ag, because Ag ^ A r . 

The interaction between the radiation standing wave, 
say 



x(x, t) = xo cos(qx) cos(wi ■ 



(21) 



and a single JV solution of the nonlocal sine-Gordon 
equation J2]), 



ip(x) = 7r + 2 arctan(Ax) 



(22) 



(both in dimcnsionless units) is well described by the 
nonrelativistic quasi-particle approach of Ref. [30l The 
JV center of mass with coordinate X(t) is subject to an 
oscillating sinusoidal trap 



w(V) + P(pVf 



pVf. 



(16) X = -[3X + xo^ e-«/ A cos(qX)cos(^ + 0), (23) 
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with an amplitude which is exponentially suppressed at 
short wavelengths, i.e., for q/X S> 1. However, for suffi- 
ciently large trap amplitudes, the vortices in each layer 
get spatially distributed with wavevector q. 

III. COMPARISON WITH NUMERICAL 
RESULTS 

The results in Refs. 0@ can be easily analyzed within 
the above theoretical framework. To make contact with 
their numerical data, one must express: all lengths in 
units of A, with A = 200 fan; the velocities in units of 
the light speed in the dielectric c = co/y 7 ?^, where cq is 
the speed of light in vacuo, and e c = 10 is the simulated 
dielectric constant; the forces in units of J/J c , where J c 
is the critical JJ current and J is the superconducting 
current across the IJJs; and the angular frequencies in 
units of the plasma gap frequency v p = uj p /2tt — c/A = 
0.47 x 10 12 Hz. Moreover, the actual layer JV density 
is p = n/L ~ 0.6 /jm" 1 with L — 100 /xm, the layer 
thickness is s — 15 A, and the penetration length ratio 
7 = 500. For this choice of numerical parameters, the 
length scales we introduced in the previous section read, 
respectively, Ag = 0.38 /im, Xem — 8.7/zm and 1/p 
0.27 fan. 

First, we note that the simulations of Refs. 00 corre- 
spond to the physical condition where Ag ^5 A r . As for 
the resonant modes xo <x e, see Eq. (fl4")l . the amplitude 
of the driving force in Eq. (f23|) turns out to scale like 
^(A/Ag) 1 / 2 , which is strong enough to drag a JV against 
the disorder field (Q} and the array of restoring forces. 
This explains the disordered spatial distribution of the 
emitting JVs, which, far from forming any ordered lat- 
tice, seem rather to get trapped by the plasma standing 
waves. In spite of the coherent nature of the plasma 
radiation, the vortex distributions in each IJJ can dif- 
fer from one another because of the intrinsic disorder 
brought about by the layer-layer coupling. 

In Fig. Q]we compare the current-voltage characteris- 
tics from simulation, reported in Fig. 4 of Ref. 0, with 
the force- velocity (f-V) curve of Eq. (j2"0)l . In the units 
of Ref. i 



- = 2.p^ P X-)=2nP v (24) 

where V = pX (V/c) is the flux- flow voltage across a 
IJJ layer and V p = v p <&q, with $o = h/2e denoting the 
flux quantum. The agreement is quite good in the linear 
regime, whereas the depinning threshold (fl9l) is clearly 
visible for J/J c — 0.2, which, in our units, corresponds 
to setting e = 1. 

The resonant plasma radiation is investigated in 
Refs. 00 on the linear branch of the J-V characteristics. 
Three plots of the plasma standing waves, two in Ref. 6 
and one in Ref.0, are shown for different J / J c ; from there 
we read out the corresponding resonance wavelengths A r . 




FIG. 1: (Color online) Current-voltage characteristics (in di- 
mensionless units) for a damping constant ft = 0.02. The 
square symbols are the simulation data extracted from Fig. 4 
of Ref. q. The red straight line is our prediction from Eq. (|24p 
for the linear Ohmic branch. The horizontal black dotted line 
is an estimate, from the numerics, of the depinning threshold 
in Eq. (19]). 
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FIG. 2: (Color online) Resonance frequency v r versus current 
intensity J (both in dimensionless units) for two values of 
the viscous constant /3. The black dots are simulation data 
extracted from Fig. 5 of Ref. 0. The two colored straight lines 
represent our theoretical predictions based on Eq. (|28[) . and 
using two different values of the damping parameter ft. 

Furthermore, the resonance frequencies v r — 2-nui r are 
either given explicitly in the text or shown in the figure 
(see, e.g., Fig. 5 of Ref. 0): The product X r v r appears 
to define a Swihart velocity, denoted here by eg, inde- 
pendent of the simulation parameters J/J c and j3, that 
is 

— ~ 0.04. (25) 
c 

In view of our radiation mechanism Q, the ratio cgjc 



6 



can be identified with c p in Eq.©; accordingly, for A ^> 1 
and 7y ~ 1, 



c 



Xp 



(26) 



is predicted to be of the order of 0.036, which is reason- 
ably close to the result in Refs. H0, given the accuracy 
of the data available. 

The average JV speed in a resonant IJJ structure is 
proportional to the resonance frequency, that is, from 
Eq. (USD, 



V 



P 



(27) 



This equation holds for all different choices of the simu- 
lation parameters presented in Refs. |6||7j. Note that the 
measured JV speeds are relatively small, V <C c, as as- 
sumed in our nonrelativistic treatment of Eq. ^ , where 
7v — 1. Moreover, when combined with Eq. (|2"4")l . this 
equation yields the dependence of v r on the simulation 
control parameters J/J c and /3. The ensuing law 



1 J 



(28) 



ductors. We show, that due to the two-scale structure 
of Josephson vortices, they radiate THz radiation on a 
characteristic scale Xem, which is much longer than the 
Josephson vortex core size Xq ~ 7s. Among all emitted 
waves, only standing modes in the sample (working as 
a cavity) survive. These standing modes produce mod- 
ulation of the density of JVs. This, in turn, make vor- 
tices mainly radiate with wavelengths corresponding to 
standing waves. Such positive feedback can result in rel- 
atively strong radiation with well pronounced maxima in 
the spectra. All our analytical estimates are in a good 
agreement with numerical data&£. 

The experimental demonstration 3 of THz radiation in 
zero magnetic field and various failed attempts at de- 
tecting THz emission in the presence of magnetic fields 
cast serious doubts on the initial idea that moving JVs 
can radiate in this frequency domain. Indeed, the now 
prevailing interpretation is that JVs ought to be con- 
sidered as perturbing degrees of freedom, which destroy 
the layer coherence and thus cause the suppression of 
THz radiation. In this study, however, we reach the con- 
clusion that, under appropriate conditions, applied mag- 
netic fields do help amplify and tune THz emission. This 
interesting result is also consistent with the recent sys- 
tematic studies in Ref. ml 



closely matches all spectral resonance peaks reported in 
Ref. |g, as shown in Fig. [2] Note that combining Eqs. 
and (|2"5)l yields the simple /3-independent relation 



Vn 



(29) 



Finally, we notice from Eqs. (HU) and $2U$ that w max 
is proportional to e 2 /(3 and V is proportional to f//3; 
as a consequence, one would expect that on decreasing 
(3 the IJJ spectral emission band shifts to lower J/J c 
while growing in intensity, both inversely proportional to 
(3. This is exactly the dependence displayed in Fig. 6 of 
Ref. S 



IV. CONCLUSIONS 

We propose a new mechanism of coherent radiation 
from the moving Josephson vortices in layered supercon- 
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